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ABSTRACT 

I discuss the recent controversy over the semileptonic branching ratio of the 
B meson, pointing out that it is only the combination of this quantity with 
the reported charm multiplicity in B decays which is in serious conflict with 
theoretical calculations. I consider possible solutions to the problem. 

1. Introduction 

There has been considerable recent interest in the question of whether there is a 
serious conflict between the measured semileptonic branching ratio of the B meson and 
its theoretical calculation.SH While a straightforward analysis of this quantity, based 
on the operator product expansion, certainly yields a result which is in conflict with 
experiment, is this discrepancy worth worrying about? In particular, is it pointing us 
toward new physics beyond the Standard Model, or telling us something interesting 
about QCD, or neither? What is the most conservative, and hence most plausible, 
explanation of the puzzle? Here I will argue that its resolution lies most likely within 
QCD itself and is probably quite mundane, namely that the calculation of inclusive 
nonleptonic decay rates is considerably less trustworthy than the calculation of inclusive 
semileptonic decays and fails even for total energy releases as large as mj. 

I will concentrate on the theoretical foundation of the calculation, explaining the 
roles which are played by the Operator Product Expansion (OPE) and the two versions 
of the duality hypothesis, "local" and "global" . We will see that if we combine the OPE 
and the assumption of local duality with the reported semileptonic branching ratio of 
the B meson, then we are led to the expectation of excess charm production in inclusive 
B decays, as much as 1.3 charm quarks per decay. The charm data, as reported, do 
not support such an excess; if anything, they indicate a small deficit of charm. What 
is more, the most plausible and least "damaging" emendations to the OPE do not 
solve this problem, as they also would imply excess charm. The real puzzle is not the 
semileptonic branching ratio alone, but rather the combination of it with the reported 
charm multiplicity. We are led to the conclusion that, if the data are confirmed, the 
key assumption of local duality may not be justified in nonleptonic B decays. 

*Invited talk presented at the Eighteenth Johns Hopkins Workshop on Current Problems in Particle 
Theory, Florence, Italy, August 31-September 2, 1994. 
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2. Theoretical Analysis of Inclusive B Decays 

The inclusive decay rate of the B meson may be organized by the flavor quantum 
numbers of the final state, 



TOT 



c tv) + r (b -^cud')+T (bices'). (1) 



We neglect rare processes, such as those mediated by an underlying b —>■ u transition 
or penguin-induced decays. We denote by d' and s' the approximate flavor eigenstates 
(d' = dcos^i — s sin 6*i, s' = dsmOx + scos^i) which couple to u and c, respectively, 
and we ignore the effect of the strange quark mass. It is convenient to normalize the 
inclusive partial rates to the semielectronic rate, defining 

r(6 -> cud') T(b -> ccs') 

^ d ~ 3 r(fe - ce p) ' Kcs 3T(b^cev)- [ ) 

The full semileptonic width may be written in terms of the semielectronic width as 

T(b^c£v) = 3f(m T )T(b^cev) J (3) 

where the factor 3f(m T ) accounts for the three flavors of lepton, with a phase space 
suppression which takes into account the r mass. Then, since the semileptonic branch- 
ing ratio is given by Br(b — > civ) = T(b c£u)/Ttot, we may rewrite Eq. ([!]) in the 
form 

f/ ,1- Br(b -> civ) 
R ud + R cs = f(m T ) Br{h \ cr , ) ■ (4) 



The measured partial semileptonic branching fractions ar 

Br(B^Xe9) = 10.7 ±0.5%, 
Br(B_^> XnD) = 10.3 ±0.5%, 

Br{B^XTv) = 2.8 ±0.6%, (5) 

leading to a total semileptonic branching fraction Br(b — ► civ) of 23.8 ± 0.9%, with 
the experimental errors added in quadrature. Of the semileptonic rate, 11% comes 
from decays to r, corresponding to a phase space suppression factor f(m T ) = 0.74, 
consistent with what one would expect in free quark decay.ii If we substitute the 
measured branching fractions into the right-hand side of Eq. , we find 

Rud + R cs = 2.37 ±0.12. (6) 

This is the result which we must try to reproduce theoretically. 

To compute R u d and R cs , it is necessary to compute individually the partial widths 
T(b — ► cev), T(b — > cud') and F(b — > ccs'). Because the computation is conceptually 
more straightforward for this case, we begin with the nonleptonic widths. 
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2.1. Nonleptonic decays 

The nonleptonic quark decay b — > cq\q 2 is mediated by four-quark operators of 
the form 

= cT 1 bq 2 T 2 q 1 , (7) 

where the denote generic Dirac matrices. Such operators are induced by the exchange 
of W bosons which have been integrated out the theory at the scale \i = Mw, where 
r"i = T 2 = 7^(1 — 75) and the coefficient of O in the effective Hamiltonian is Gf/V^- As 
the theory is evolved down to the scale /1 = m&, this coefficient changes, and additional 
operators of the form (|7|) are induced with coefficients which depend on Xog^/m^). In 
this schematic discussion of the inclusive decay rate, we will ignore the interference 
between these operators and treat only the generic operator O above. 

The inclusive decay rate of the B meson mediated by the operator O is given by 

T ^^(B\0^\X)(X\0\B) = J2\(X\0\B)\\ (8) 

x x 

where the sum runs over all states X with the appropriate quantum numbers. This 
expression may be rewritten via the optical theorem in terms of the imaginary part of 
a forward scattering amplitude, 

T ~ Im (B\ T{0\ O} \B) . (9) 

The next step is to compute the time-ordered product T{0\ O} in perturbative QCD. 

Before we proceed with this step, it is useful to examine the conditions under 
which it is sensible to apply a perturbative analysis. In Fig. [I], I present a cartoon of 
the quantity 

f(m Q )=Im(M Q \T{0\0}\M Q ), (10) 

where the b quark has been replaced with a quark Q of mass rriQ and the B meson 
by the ground state meson Mq. The expectation value of the correlator is given as a 
function of mq. We see in Fig. p] that the behavior of /(ttiq) is qualitatively different 
for rriQ large and for rriQ small compared to the scale Aqcd which characterizes the 
nonperturbative regime of QCD. For small tuq, the correlator is dominated by a few 
resonances whose position and strength depends on nonperturbative strong dynamics. 
In this "resonance regime" , no perturbative calculation can reproduce the rapid varia- 
tions in the shape of /(mg). For large rriQ, however, these rapid variations damp out, 
as in this region the correlator is dominated by multiparticle states. While even here 
there are new thresholds associated with the production of additional pions, their effect 
is small compared to the smooth background of states to which they are being added. 
In this regime, one might expect the smooth output of a perturbative calculation of the 
correlator to approach the physical fimq). This property, that at large enough energies 
one may compute accurately within perturbative QCD, is known as "local duality". 
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Fig. 1. Cartoon of f(m Q ) = Im (M Q \ T{0\ O} \M Q ). 



Where is the transition between the resonance and perturbative regimes, or in 
other words, above what value of itlq does one expect local duality to hold? In particu- 
lar, is mt, sufficiently large for perturbative QCD to be applicable? Since we cannot vary 
nib experimentally, we cannot answer this question by tracing out f(rriQ) and looking 
for where it becomes smooth. We can only argue on dimensional grounds. Naively, we 
might expect to be well in the perturbative regime, since rrib fa 5 GeV is large compared 
to typical QCD scales of hundreds of MeV. However, the situation is not necessarily 
so clear, since the nonleptonic decay is into three colored particles which must divide 
the available energy between them, and charm quark(s) in the final state carry away 
substantial rest energy as well. We also note that the transition point between the 
resonance and perturbative regimes undoubtedly depends on the operator O, as well 
as on the external states Mq. 

While we are free to have our prejudices, it is probably best to leave this question 
open for the time being. We may plunge ahead with the perturbative calculation of 
the total width, while bearing in mind that the assumption of local duality is made by 
necessity but cannot rigorously be justified. Later, it will be necessary to evaluate any 
discrepancies between experiment and the theoretical results in light of their sensitivity 
to this crucial assumption. 

It is also true that local duality may well fail at different p 2 in the channels 
b — > cud' and b — > ccs'. In fact, we might expect it to fail first (at higher energy) in the 
channel with two charm quarks in the final state, since, due to the charm quark rest 
energy, the kinetic energy released per particle is lower in this channel. For example, 
in a simple free quark decay model the s quark has an average energy of only 1 GeV. 
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This intuitive expectation is our first hint, with more to come, that the quantity R u d 
is more reliably calculable than R cs . 

2.2. Semileptonic decays 

The situation is considerably better for semileptonic decays. Here the decay is 
mediated by a product of currents, 

= cr(l-j 5 )b£ lfM (l-j 5 )u, (11) 

and the matrix element factorizes into hadronic and leptonic pieces, 

(X £( Pe )u(p,)\ J£j ltl \B) = (X\ r h \B) ■ (e(p £ )9(p,)\J^\0) . (12) 

We then find an expression in which the integral over the momenta of the leptons is 
explicit, 

r ~ J dy dv ■ q dq 2 L^v ■ q, q 2 , y) W^(v ■ q, q 2 ) , (13) 

where L^ v is the lepton tensor and W^ u the hadron tensor. Here the momentum of the 
external h quark is written as p b = . The other independent kinematic variables 
are q^ = p^ + pp and y — 2v ■ pe/m b . It is convenient to scale all momenta by m b , so 
q = q/rrib. The hadronic tensor is given by 

w>" = J2(B\j?\x)(x\r h \B) 

x 

= -21m (B\i J dxe^T{4\x),r h (0)} \B) 

= -2ImT^. (14) 

We may perform the integrals in y, v ■ q and q 2 in Eq. ([13|) to compute the total 
semileptonic decay rate, or leave some of them unintegrated to obtain various differen- 
tial distributions. 

Let us consider as an example the doubly differential distribution dT/dydq 2 , for 
which we must perform the integral over v-q, for y and q 2 fixed. The range of integration 
is given by (y + q 2 / y)/2 < v ■ q < (1 + q 2 — m 2 )/2, where m q is the mass of the quark to 
which the b decays, and m q = m q /m b . The integration is pictured as the solid contour 
in Fig. ||, along with the analytic structure of T M ^; the cut extends along the real axis 
except in the region (1 + q 2 — rh 2 )/2 < v ■ q < {{2 + m q ) 2 — q 2 — l)/2. Note that we have 
already included only the imaginary part of T^ v by integrating over the top of the cut 
and then back underneath. In general, along the cut T^ u will look qualitatively like 
the cartoon shown in Fig. [IJ depending on v ■ q in a complicated nonperturbative way. 
However, we may now use Cauchy's Theorem to deform the contour of integration to 
the dashed contour in Fig. [|, which lies far away from the cut everywhere but at its 
endpoints. Along this new contour, we are far from any physical intermediate states, 
and it is safe to perform the operator product expansion in perturbative QCD. 
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Fig. 2. Cuts and contours in the complex v ■ q plane, for fixed q 2 and y. The cut extends 
along the real axis except in the region (1 + q 2 — m 2 ) /2 < v • q < ((2 + fn q ) 2 — q 2 — l)/2. The 
endpoints of the contour integral are at v ■ q = (y + q 2 /y)/2 ± ie. 

The additional integrals over the lepton phase space allow us to determine the 
physical rate in terms of certain smooth integrals of T^ u rather than in terms of the 
time-ordered product evaluated at a single point. These smooth integrals may be com- 
puted reliably by deforming the contour into the unphysical region. That we we can 
compute integrals of T^ v in this way is the property of "global duality" . In principle, 
it rests on a much firmer theoretical foundation than does local duality, which asserts 
the computability of QCD correlators directly on the physical cut rather than off in 
some unphysical region of the complex plane. 

Unfortunately, the dashed contour must still approach the physical cut near the 
endpoints of the integration, introducing an unavoidable uncertainty into the calcula- 
tion. There is a two-pronged argument that this uncertainty is small. First, for large 
rrib, the proportion of the contour which is within Aqcd of the physical cut scales as 
and thus makes a small contribution to the total integral. Second, in this 
small region one may invoke local duality, as in the case of nonleptonic decays. The 
endpoints of the integration lie conveniently at the point of maximum recoil of the final 
state quark, where this property is most likely to hold. In addition, since there are fewer 
colored particles in the final state than in nonleptonic decays, one would expect this 
assumption to be more reasonable in the semileptonic than in the nonleptonic case. 

The point is that semileptonic decays, because they rely on global rather than 
local duality, may be calculated much more reliably in perturbative QCD than may 
nonleptonic decays. This is why we have chosen, in defining R u d and R cs , to normalize 
the nonleptonic widths to the semielectronic width. (Of course, this also eliminates 
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Fig. 3. The operator product expansion of T{0^,0}. 

the dependence on the CKM parameter V c b and some of the dependence on m&.) The 
theoretical uncertainties in R u d and R cs will be dominated entirely by the uncertainties 
in the numerators. 

2.3. The perturbative calculation 

The computation of TUDj. O} in perturbation theory follows well-known tech- 
niques for both semileptonic&S and nonleptonic0 decays. Since the quarks in the 
intermediate state are assumed to be far from the mass shell by an amount of order 
m.&, we may develop the operator product in a simultaneous expansion in l/m& and a s . 
Graphically, the expansion proceeds as in Fig. |3] in the nonleptonic case, to whatever 
accuracy is desired. The Operator Product Expansion takes the form of an infinite 
series of local operators, 

T{0\0} = Y,^Om (15) 

n=0 m b 

where the operator O n has dimension n + 3 and the c n are calculable functions of a s . 
The notation is somewhat schematic, in that more than one operator may appear for a 
given n. For example, the operators of lowest dimension which appear in the expansion 
are of the form 

C = bTb, 

2 = bT^D^b, 

where D„ is the covariant derivative, and r^ " are combinations of Dirac matrices. 



7 



Once the OPE has been performed, we must compute the expectation values of 
the local operators O n in the B state. Such expectation values involve the structure 
of the meson as determined by nonperturbative QCD, and typically are incalculable. 
However, here we may use the additional symmetries of the heavy quark limitJl! to 
simplify the description of these matrix elements. In this limit, we construct an ef- 
fective theory in which the production of heavy quark-antiquark pairs is eliminated, 
at the price of introducing into the Lagrangian an infinite tower of nonrenormalizable 
operators suppressed by powers of l/m&. In the effective theory, we may compute the 
matrix elements which arise at a given order in n in terms of a few nonperturbative 
parameters. In fact, the matrix element of Oq is determined exactly by the conservation 
of heavy quark number; if the expansion is truncated at this point, then we reproduce 
precisely the result of the free quark decay model. Since all matrix elements of 0\ van- 
ish,! the leading corrections to free quark decay come only at relative order l/mf, not 
at order l/m?,. These corrections come from initial state binding effects of the b quark 
in the B meson, the energy of the heavy quark being altered from its "free" value by 
its interactions with the cloud of light quarks and gluons which surrounds it. They are 
governed by only two parameters, defined by the expectation values 



The parameter Ag, associated with an operator which violates the heavy quark spin 
symmetry ,0 is determined by the mass difference between the B and B* mesons, 



Unfortunately, Ai cannot at this time be measured experimentally, and for its esti- 
mation we must rely on models of QCD. The same is true of the matrix elements of 
operators of dimension six and higher. 

The ratios R u d and R cs receive contributions, then, both from perturbative (a s ) 
and nonperturbative (Ai, A2, • • .)^ources. I will summarize below the analysis of the 
semileptonic"^ and nonleptonid^HJ radiative corrections which exists in the litera- 
ture. The nonperturbative correct ions&Hl! to R u d and R cs have recently been studied 
up to operators of dimension six and shown to be quite smallS; because they turn out 
to be unimportant compared to the uncertainties in the radiative corrections, I will not 
include them from here on. 

After some manipulation,! the ratios R u d and R cs may be written in the form 



(B\ b{iD) 2 b\B) =2M B X 1 , 

(B\ b{-\a^)G^ v \B)=QM B X 2 . 



(16) 




(17) 



R 



'CS 



P(/i) + SP ud (/i, m c ) , 
G(m c /m b ) P(/i) + <5P cs (/i, m c ) 



(18) 



Here 



P(/i) = ?7(/i) 1 + 



a s (n) 



(19) 



7T 



8 



1 2 3 4 5 

(5) 

Fig. 4. The radiative correction P{[i). The upper curve corresponds to A^ = 220 MeV, 
the middle curve to A^- = 180 MeV, and the lower curve to A^- = 140 MeV. We take 

MS ' MS 

m b = 4.8 GeV. 

is the radiative correction in the massless (m c = 0) limit, with leading logarithmic!!! 
(r]([i)), one loop0 (a s (/i)/7r) and subleading logarithmic^ (J 2 (/i)) contributions. The 
inclusion of the subleading logarithms softens the dependence of -P(yu) on the renor- 
malization scale, which is shown in Fig. |] for a range of values of A|^. The quantities 
SP u d and 5P CS parameterize the effect on the radiative corrections of the charm quark 
mass. Recent calculations^ show 5P u d to be no larger than ~ 0.02. There are 
indications&0 that 5P CS may be substantially larger than 5P u d, but a full calculation 
is not yet available; we consider the uncertainty due to SP CS below. 

While there is certainly some variation in R u d and R cs due to the dependence of 
the radiative corrections on the renormalization scale, the largest uncertainty in R cs 
comes from the phase space factor G{m c jm^). There is an additional dependence on the 
charm mass in this channel because of the second charm quark in the final state. The 
analytic expression for G{m c /m b ) has been computed!; we present a plot in Fig. ||. Here 
we have used the fact that within the heavy quark expansion, the difference between 
m c and m;, is much more precisely known than either of the masses individually, and 
is given in terms of the spin-averaged D meson and B meson masses: 

m b -m c = (M B ) avc . - (M D ) avc . = 3.34 GeV . (20) 

In Fig. |5| we hold m b — m c fixed and consider variations of mj only. 

The important inputs into the theoretical calculation, then, are /x, m b and, to a 
lesser extent, A^. As first emphasized by Bigi et al.,0 it is not a trivial task to fix 
these parameters to satisfy the experimental constraint R u d + R cs = 2.37 ± 0.12. For 
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Fig. 5. The phase space suppression factor G(m c /mb), as a function of m;, with m;, — m c = 
3.34 GeV held fixed. 

(5) 

example, if one takes the reasonable values /i = m;, = 4.8 GeV and = 180 MeV 
and neglects 5P u d and 5P CS , then -P(/i) = 1.27, G{ra c jm^) = 0.36 and R u d + R cs is only 
1.73, far short of the experimental number. The question, then, is whether there exist 
values of the inputs for which we can manage to reproduce the data. 

It is useful to start by considering R ud , for which the calculation is likely to be 
more reliable. To enhance this quantity we may vary the renormalization scale fi. The 
choice /i = m b is motivated by the fact that the total energy released in the decay b — > 
cud' is mb] however, since this energy has to be divided among three colored particles, 
perhaps the appropriate scale is lower.Q For fi = 1.6 GeV ~ rrib/3, for example, and 
A^|t = 220 MeV, we have R u d = -P(/i) = 1.52, a modest enhancement. If we take this 
value of /i as a reasonable lower limit, then the data on the semileptonic branching 
ratio implies R cs > 0.85. This restriction on R cs follows simply from the assumption 
that the perturbative calculation of R u d is reliable. 

There is more than one way to achieve R cs > 0.85, both within the present calcu- 
lation and by going outside it. The easiest, perhaps, is to lower irtb to 4.4 GeV, which 
corresponds to m c = 1.1 GeV. These masses for the bottom and charm quarks are some- 
what low, but probably not unacceptably so. Taking \i = 1.6 GeV, then, we would find 
G(m c /m b ) = 0.58, and R cs = 0.89. Alternatively, we note that recent calculations&0 
indicate that 5P CS may not be negligible at all, and that the inclusion of charm mass 
effects in the radiative correction in this channel may lead to an enhancement of as 
much as thirty percent .0 Finally, we could postulate the failure of local duality in the 

*This scenario is_£.upported by a recent estimateE] of /i for semileptonic decays using the BLM scale- 
setting criterion£3 in which a scale fj, as much as an order of magnitude below mi, is indicated. 
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b — > ccs' channel, since this is the channel in which we expect this assumption to be 
the least reliable. From any of these points of view, or from a combination of them, 
the B semileptonic branching ratio is surprising, perhaps, but certainly explicable. We 
simply have to allow for an enhancement of the final state with two charmed quarks 
relative to the final state with one, as compared to our naive expectation. At the same 
time, such an enhancement is absolutely necessary if the data are to be reconciled with 
the theoretical prediction for R ud . 

3. Experimental consequences 

This proposal that an enhancement of R cs explains the observed semileptonic 
branching ratio is not without experimental implications of its own. The most striking 
of these is a large number of charmed quarks per B decay. The charm multiplicity n c 
may be written 

BrfB^XJu) , . 

n c = 1 + R cs \, , ; . (21) 

Using Br(B -> XJv) = 23.8% and f(m T ) = 0.74 in Eq. (0) yields 

n c = 1.00 + 0.32 R cs , (22) 

which for the values of R cs necessary to explain the semileptonic branching ratio would 
indicate n c ~ 1.3. 

While there has been no direct inclusive measurement of n c , it may be estimated 
by summing individual partial widths. Among the B decay products, there are contri- 
butions to n c from charmed mesons, charmed baryons, and cc resonances. The number 
of charged and neutral D mesons per decay, summed over B and B, has been measured 
to be 



n D ± = 0.246 ±0.031 ±0.025, 
n D0 - = 0.567 ±0.040 ±0.023. (23) 

The branching ratio to Df mesons has not yet been determined, because no absolute 
D s branching ratio has been measured. However, it is known 

n D ± = (0.1181 ± 0.0053 ± 0.0094) 

and the branching ratio for D s — > (jnt is expected to be about 3.7%. Inclusive B decays 
to cc resonances below DD threshold contribute two units to n c . The branching ratio 
to ip has been measured to be (1.11 ± 0.08)%, including feed-down from the excited 
states ip' an d Xc-^ Combined with the measurementsi Br(B ip'X) = (0.32±0.05)%, 
Br[B -> XdX) = (0.66±0.20)% and Br[B -> r} c X) < 1%, we expect that the inclusive 
B branching ratio to charmonium states below DD threshold is between two and three 
percent. 



3.7% 



Br(D s -> 0tt) 



(24) 
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Finally, there is the contribution to n c of B decays to charmed baryons. While the 
production of charmed baryons may be estimated from the measured inclusive produc- 



the underlying production process and the predominant decay chains.^ Such assump- 
tions introduce considerable uncertainty. Recently, however, the CLEO Collaboration 
has reported certain exclusive branching ratios of B mesons to charmed baryonsB 



with a branching ratio to A c roughly twice that to S c . Only statistical errors are included 
here. (The new CLEO data do not support the recent suggestion^! that charmed baryon 
production in B decays is dominated by the underlying process b — > ccs', which would 
yield final states with two charmed baryons rather than one, somewhat alleviating the 
"charm deficit".) 

Summing over all the contributions, then, the experimental charm multiplicity is 
n c w 1.1 ± 0.1, where the error is dominated by the uncertainty on charmed baryon 
production. Such a value does not provide evidence for the excess charm production 
indicated by the data on the semileptonic branching ratio. 

4. Theoretical implications 

What are we to make of this situation? We have seen that given the uncertainties 
in the calculation, neither the semileptonic branching ratio of the B meson nor the 
charm multiplicity in B decays is by itself much of a problem. However, the combination 
of the two measurements is in sharp contradiction with current theoretical calculations. 
While there is considerable room in the computation of either quantity to make it agree 
with experiment, the adjustments that make one quantity better make the other worse. 
When viewed in this light, there are important implications for the theory of inclusive 
B decays. 

Of course, it is certainly possible that in the future the data may move in the 
direction of the theoretical calculations, and that the discrepancy simply will go away. 
In particular, this may happen once the contributions to n c of charmed baryons is better 
understood. If future measurements give support to the prediction n c ~ 1.3, then the 
resolution of the semileptonic branching ratio puzzle simply will be that R cs is larger 
than expected. Whether this is due to a failure of local duality, to a low bottom quark 
mass, or to a renormalization scale much lower than the naive choice fi = nib, may 
still be unclear, but as we have stressed, such an enhancement is not hard to arrange 
theoretically. 

If, on the other hand, the present measurement of n c is confirmed, then the 
consequences are more extreme. The most natural explanation of such a situation would 
be that, for some reason, the calculation of R u d is also not reliable. Again, this could be 
due simply to a failure of local duality: perhaps mj is not large enough to calculate any 
inclusive nonleptonic decay widths at all. Such a conclusion would be disappointing, 




1.5 ±0.7%, 
2.4 ±1.3%, 



(25) 
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of course, since this is the nonleptonic channel in which local duality was considered 
more likely to hold. However, disappointing is not the same as unnatural, and it is 
worthwhile to recall that the invocation of local duality was based essentially only on 
"educated hope" . There is nothing profound, necessarily, in learning that this time our 
luck was not so good. 

An alternative explanation would be that the initial decomposition ([]]) of the 
total width was incomplete, in that charmless final states were neglected. A charmless 
branching ratio of 20% would be sufficient to solve the semileptonic branching ratio 
puzzle without inducing an unacceptable charm multiplicity. Of course, the neglect 
of charmless final states is well justified within the Standard Model, in which b — *■ u 
weak transitions are suppressed relative to b — ► c by a factor \V u b/V c b\ 2 ~ 0.01, and 
transitions such as b — > sg and b — > sj are found at the 10~ 3 level or below. However, 
in extensions of the Standard Model this need not be the case, and it is possible 
to construct scenarios@ in which, for example, the rate for b — > sg is substantially 
enhanced without inducing dangerous new contributions to b — > 57, for which there are 
already stringent limits. Models with an enhancement of 6 — > sg are also constrained by 
the known limits on the exclusive decay B — > Kir, although unknown strong interaction 
effects preclude one from translating these limits into a firm bound on the strength of 
the b — > sg coupling. 

Although it is certainly exciting to consider such novel possibilities, I must stress 
that they are not as yet well motivated by the data. The most conservative, natural 
and plausible resolution of the semileptonic branching ratio puzzle lies within QCD 
itself. In fact, the entire subject ought best to be thought of as an interesting test of 
the applicability of the Operator Product Expansion and perturbative techniques to 
the calculation of inclusive nonleptonic decays in the B system. As of today, the jury 
remains out on this important question. 
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